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ABSTRACT 

We compare X-ray hydrostatic and wcak-lensing mass estimates for a sample of 12 clusters that 
have been observed with both XMM-Newton and Subaru. At an over-density of A = 500, we obtain 
1 — M X /M WL = 0.01 ± 0.07 for the whole sample. We also divided the sample into undisturbed and 
disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation 
parameters, obtaining 1 - M X /M WL = 0.09±0.06 and -0.06±0.12 for the undisturbed and disturbed 
clusters, respectively In addition to non-thermal pressure support, there may be a competing effect 
associated with adiabatic compression and/or shock heating which leads to overestimate of X-ray 
hydrostatic masses for disturbed clusters, for example, in the famous merging cluster A1914. Despite 
the modest statistical significance of the mass discrepancy, on average, in the undisturbed clusters, we 
detect a clear trend of improving agreement between M x and M WL as a function of increasing over- 
density, M X /M WL = (0.908 ± 0.004) + (0.187 ± 0.010) ■ log 10 (A/500). We also examine the gas mass 
fractions, / gas = M gas /M wh , finding that they are an increasing function of cluster radius, with no 
dependence on dynamical state, in agreement with predictions from numerical simulations. Overall, 
our results demonstrate that XMM-Newton and Subaru are a powerful combination for calibrating 
systematic uncertainties in cluster mass measurements. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: clusters: indi- 
vidual (Abell 1914) — gravitational lensing: weak — surveys — X-rays: galaxies: 
clusters 
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The mass function of galaxy clusters depends on both 
the matter density and the expansion history of the uni- 
verse. Indeed, clusters provided early evidence for a low 
density universe (White et al. 1993). Clusters, as pow- 
erful tools to constrain cosmological parameters (e.g., 
Zhang & Wu 2003; Balogh et al. 2006; Henry et al. 

2009) , currently receive much attention as a potential 
probe of the dark energy equation of state parameter 
(w = p/p, where p is the energy density and p is the pres- 
sure), through the evolution of the mass function, (e.g., 
Vikhlinin et al. 2009a, 2009b). In addition, gas mass 
fraction measurements also potentially provide an impor- 
tant cosmological probe, under the assumption that gas 
mass fractions do not evolve with redshift (e.g., Vikhlinin 
et al. 2002; Allen et al. 2004, 2008; Mantz et al. 

2010) . Upcoming galaxy cluster surveys will shortly de- 
liver huge amounts of multi- wavelength data, e.g., from 
Subaru/Hyper-Suprime-Cam, eROSITA, PLANCK, and 
South Pole Telescope (SPT). To achieve good control 
over systematic errors in cosmological measurements, for 
example of w, based on these surveys, it is crucial to 
understand cluster mass estimates. 

The density and temperature distributions of the hot 
X-ray emitting gas within galaxy clusters can be used 
to estimate the total mass of the cluster. Since the ac- 
celeration of the gas is still not well understood, it is 
assumed that the acceleration terms are negligible when 
estimating cluster masses from X-ray data; the resulting 
mass estimates thus invoke the assumption of hydrostatic 
equilibrium (H.E.), and are hereafter referred to as "X- 
ray hydrostatic mass estimates" . Such mass estimates 
also assume that the total pressure is dominated by the 
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thermal pressure of the gas. 

Numerical simulations (e.g., Evrard 1990; Lewis et al. 
2000; Rasia et al. 2006; Nagai et al. 2007; Piffaretti & 
Valdarnini 2008, PV08, hereafter; Jeltema et al. 2008; 
Lau et al. 2009) have pointed out that X-ray hydrostatic 
mass estimates may underestimate cluster mass, and that 
this effect is most pronounced for clusters classified as 
disturbed, based on their X-ray morphology. There- 
fore, measurements of cosmological parameters based on 
the X-ray-measured cluster mass function or the X-ray- 
measured gas mass fraction may be biased. It is thus 
crucial, as a minimum, to calibrate observationally both 
the X-ray-measured cluster masses and gas mass frac- 
tions for a representative cluster sample. This opportu- 
nity is offered by measuring cluster masses using both 
X-ray and gravitational lensing data. 

X-ray and gravitational lensing mass measurements 
have complementary advantages and disadvantages. The 
n\ dependence of the X-ray emissivity of the intraclus- 
ter medium (ICM) helps guard against projection effects 
due to mass along the line of sight through the cluster. 
However, as discussed above, X-ray analysis requires as- 
sumptions to relate electromagnetic radiation to a model 
of the mass distribution. An ideal case is to use an ap- 
proach that is insensitive to the dynamical state to esti- 
mate the cluster mass. Gravitational lensing fulfills this 
requirement because the lensing signal is insensitive to 
the physical state and nature of the deflecting matter 
distribution. Lensing is, however, prone to projection 
effects because it is sensitive to all mass along the line 
of sight through the cluster (e.g., Hoekstra 2001). It 
is thus of paramount importance to combine these two 
techniques to develop a thorough understanding of clus- 
ter mass measurements. Since the 1990's, the use of X- 
ray and lensing observations has been proposed to test 
deviations from H.E. and extra pressure support (e.g., 
Miralda-Escude & Babul 1995; Wu & Fang 1996; Squires 
et al 1996; Allen 1998; Zhang et al. 2005, 2008; Mahdavi 
et al. 2008; Richard et al. 2010). Most earlier studies 
employed "target of opportunity" mode — i.e., consid- 
ered clusters with available data — making the studies 
biased by design. 

The Local Cluster Substructure Survey (LoCuSS 17 ; G. 
P. Smith et al. in preparation; Zhang et al. 2008; Haines 
et al. 2009a, 2009b; Sanderson et al. 2009; Marrone 
et al. 2009; Okabe et al. 2010a; Richard et al. 2010) 
is a systematic multi-wavelength survey of X-ray lumi- 
nous (£x,o.i-2.4kov>2 • 10 44 erg s _1 ) galaxy clusters at 
0.15 < z < 0.3 selected from the ROSAT All-Sky Survey 
(RASS; Ebeling et al. 2000; Bohringer et al. 2004) in a 
manner blind to the dynamical status. As a first step to- 
ward a comprehensive X-ray /lensing study we assembled 
a sample of 19 clusters with archival XMM-Newton ob- 
servations, and weak-lensing mass measurements in the 
literature (Zhang et al. 2008). The mean weak-lensing 
mass to X-ray hydrostatic mass ratio was found to be 
(M WL /Af x ) = 1.09 ± 0.08 at an over-density of A = 500 
with respect to the critical density. Mahdavi et al. (2008) 
found the average X-ray to weak-lensing mass ratio to 
be 0.78 ± 0.09 at A = 500, being consistent with the 
results in Zhang et al. (2008). The uncertainties in the 
above analysis was dominated by measurement errors, in 

17 http: / /www. sr. bham.ac.uk/locuss 



particular, the lensing masses drawn from the literature 
based on early lensing data (Bardeau et al. 2007; Dahle 
2006), with typical seeing of FWHM^0."8, using four dif- 
ferent cameras on three different <3.6-m telescopes, with 
fields of view (FOVs) spanning r~0.7— 2 Mpc at z ~ 0.2. 
The next step in our program is to employ uniform high- 
quality weak-lensing data from our own dedicated ob- 
serving program with Subaru/Suprime-CAM. Recently, 
the first batch of weak-lensing mass measurements have 
become available for 22 clusters based on data taken in 
good conditions (FWHM~0."6) through two-filters with 
details outlined in Section l2~T1 (Okabe et al. 2010a; Okabe 
<fc Umetsu 2008). These clusters were selected based on 
observability from Mauna Kea on nights allocated to Lo- 
CuSS, and are thus unbiased with respect to their X-ray 
properties. Archival XMM-Newton data are available for 
12 of the 22 clusters — see Zhang et al. (2008) for de- 
tails. In this work, we therefore compare X-ray (based on 
the assumption of H.E.) and weak-lensing mass estimates 
for these 12 clusters, and also compare our observational 
results with predictions from numerical simulations. It 
is worth noting that Subaru/Suprime-CAM provides an 
excellent match to the XMM-Newton FOV, and also that 
these Subaru weak-lensing mass estimates have optimal 
precision in the density contrast range of 500 < A < 2000 
(Okabe et al. 2010a), again, well matched to the XMM- 
Newton X-ray estimates. 

The outline of this paper is as follows. In Section [2] 
we briefly describe the weak-lensing and X-ray analysis. 
The results are presented in Section |3l and discussed in 
Section 2J We summarize our conclusions in Section [SJ 
Throughout the paper, we assume O m = 0.3, il\ = 0.7, 
and Hq = 70 kms _1 Mpc _1 . Confidence intervals cor- 
respond to the 68% confidence level. Unless explicitly 
stated otherwise, we apply the Orthogonal Distance Re- 
gression package (ODRPACK 2.01 18 , e.g., Boggs et al. 
1987) taking into account measurement errors on both 
variables to determine the parameters and their errors of 
the fitting. 

2. DATA ANALYSIS 
The 12 clusters in our sample are listed in Table [1] 

2.1. weak-lensing analysis 

The Subaru %'- and l/-band imaging data and detailed 
weak-lensing analysis are described in Section 2 in Ok- 
abe et al. (2010a), in which the i'-band data are used 
for shape measurements, and the U-band data are used 
to remove foreground and cluster galaxies. Excluding 
unlensed galaxies from the background galaxy catalog is 
of prime importance for accurate weak-lensing mass esti- 
mates, especially at higher over-densities. This so called 
dilution bias increases as a function of the density con- 
trast A, and can cause and M^Q^ to be biased 
low by - 20% - 50% (Okabe et al. 2010a). Our two- 
filter lensing data allow us to construct secure red+blue 
background galaxy samples, defined as faint galaxies with 
colors that are redder and bluer than the cluster red- 
sequence by a minimum color offset. This strategy typi- 
cally reduces the dilution bias to per cent level (Okabe et 
al. 2010a). The mean redshift for the background galaxy 

18 http://www.netlib.org/odrpack and references therein 
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catalog is estimated by matching the magnitudes and col- 
ors of background galaxies to the COSMOS photometric 
redshift catalog (Ilbert et al. 2009). More precisely, the 
mean redshift is computed as a lens weighted average 
over the redshift distribution, dPy/h/dz, which is defined 
by (D ls /D s ) = J Zd dzdPwh/dzDi s /D s , where D s and As 
are the angular diameter distances to source and between 
lens and source, respectively. 

Galaxy cluster mass distributions are often modeled 
as Navarro, Frenk, and White (NFW; Navarro ct al. 
1997) halos or singular isothermal spheres (SIS). In weak- 
lensing studies, the distortion signal is used to constrain 
the model parameters. Okabe et al. (2010a) have shown 
that the NFW model fits the lensing distortion profiles 
well in both statistical studies and their analysis of in- 
dividual clusters. The SIS model is statistically inade- 
quate to describe stacked tangential distortion profiles 
with pronounced radial curvatures, and is rejected at 
6(7 and lltr level, respectively, for their two sub-samples 
with the virial masses in the ranges of < 6 x 1O 14 /i _1 M0 
and > 6 x 10 14 h~ 1 MQ. It is also important to note that 
the mean ratio of masses obtained from SIS and NFW 
models is 0.70±0.05 at r^Q for our sample of 12 clusters. 

The projected mass distribution can also be obtained 
using a model-independent approach, the so-called £ c - 
statistics method (Fahlman et al. 1994; Clowe et al. 
2000), which is complementary to the tangential shear fit 
method. The £ c -statistics method measures the discrete 
integration of averaged tangential distortions of source 
galaxies outside given radii. The £ c -statistics method is 
thus less sensitive to the detailed structure of clusters on 
small scales than using models fitted to the full tangential 
shear data. Okabe et al.'s Cc-based model-independent 
projected masses are in good agreement with the pro- 
jected masses computed from the NFW-based models at 
A = 500, but not with the SIS-based models. 

We also compared the NFW-based spherical mass es- 
timates for our sample of 12 clusters with the spheri- 
cal mass estimates based on deprojecting the Cc-based 
model-independent masses. In the latter case, the spher- 
ical mass was obtained by assuming an NFW profile. The 
spherical masses using the (^-statistics and NFW mod- 
els are statistically consistent, given their mean ratios of 
1.00 ± 0.08, 0.96 ± 0.07, and 1.02 ± 0.07, at r^o , r^ , 
and r^Q Q , respectively, weighted by the inverse square 
of the errors. For two clusters, A383 and A2390, the 
masses from the NFW models are lower than the de- 
projected masses from the £ c -statistics method due to 
substructures in the cluster cores. 

Based on the above tests, our joint analysis uses the 
NFW masses obtained from the tangential distortion pro- 
files in Okabe et al. (2010a). The three-dimensional 
spherical cluster masses (see Table [2]), M^ L , within a 
sphere of radius r&, is derived following the NFW model, 
p oc r (1 + ca^/^a) 2 , where ca is the concentration 
parameter. 

The weak-lensing analysis includes both statistical er- 
ror and errors of the photometric redshifts of source 
galaxies in mass measurements. The latter is estimated 
by bootstrap re-sampling to match the COSMOS cata- 
log. Leauthaud et al. (2010) found that the latter mainly 
matters for galaxies close to the lens. The typical la 
total uncertainty on the weak-lensing mass estimates is 



~ 12% — 21% at A = 500 with four exceptions, i.e., 
RXJ2129.6+0005, 39%; A68, 41%; A115 (south), 53%; 
and Z7160, 42%. 

Hoekstra (2001) pointed out that the uncertainty in 
weak-lensing mass estimates of clusters, caused by dis- 
tant large-scale structures (uncorrelated) along the line 
of sight is fairly small for deep observations (20 < R < 
26) of massive clusters at intermediate redshifts. The 
typical lcr relative uncertainty is about 6% if the lensing 
signal is measured out to 1.5h^ Mpc. All 12 clusters in 
our sample are massive clusters (> 5 keV) at intermedi- 
ate redshifts (0.15 < z < 0.3) with deep Subaru obser- 
vations to ~ 26 mag. The Subaru FOV covers the entire 
cluster up to a few Mpc for our sample. Therefore, the 
mean mass uncertainty caused by large-scale structures 
due to projection should be no more than 6%. Therefore 
we neglect this error in the mass estimates. 

2.2. X-ray analysis 

The X-ray analysis was carried out independently from 
the weak-lensing analysis. The description of XMM- 
Newton data and X-ray mass modeling is given in Sec- 
tion 2 in Zhang ct al. (2008), and more details on the 
data reduction are given in Section 2 in Zhang et al. 
(2006, 2007). The radial temperature profile was mea- 
sured using spectral data including deprojection. The 
X-ray spectrum measuring the global temperature was 
used to calculate the cooling function for the conversion 
from X-ray surface brightness profile to electron num- 
ber density profile, in which we included both depro- 
jection and XMM-Newton point-spread-function correc- 
tions. The gas mass M gas (< r^ L ) was derived by inte- 
grating the electron number density, which was fitted by 
a double-/? model, and assuming /i e = 1.17. The X-ray 
hydrostatic mass M x (< r^ h ) was measured from the 
temperature profile and electron number density profile 
assuming spherical symmetry and H.E.. Unless explicitly 
stated otherwise, both the X-ray gas mass and the hydro- 
static mass are measured within the radius r^ L obtained 
in the weak-lensing analysis. The typical la uncertainty 
on the X-ray hydrostatic mass is ~ 18%— 30% at A = 500 
with no exceptions (see Table [2]). 

2.3. Selection of cluster centers 

The cluster centers were derived independently in the 
weak-lensing and X-ray analysis. In the weak-lensing 
analysis, we follow Okabe ct al. (2010a, see their Sec- 
tion 3.2) and adopt the position of the brightest cluster 
galaxy (BCG) as the cluster center. This is motivated 
by the coincidence of the peak of the two-dimensional 
cluster mass distribution with the BCG position, and 
the strong gravitational lensed images centered close to 
the BCG position in some of the clusters (Richard et al. 
2010). 

In the X-ray analysis in Zhang et al. (2008), the clus- 
ter center is determined using the flat fielded X-ray image 
in the 0.7-2 keV band as follows. The procedure is ini- 
tiated by deriving the first flux-weighted center within 
a 1' aperture centered at the peak of the cluster X-ray 
emission. Iteratively, we re-derive the flux-weighted cen- 
ter but within the aperture, which is 1' larger than the 
previous one and centered at the previous flux-weighted 
center, till the coordinates of the flux- weighted center do 
not vary anymore. The iteration is less than 10 times to 
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fulfill the goal. The final flux-weighted center is taken as 
the X-ray cluster center. 

We list the lensing and X-ray centers in Table [TJ they 
agree to within O.f4r25oo for all except two clusters, 
namely A1914 and RXCJ2337.6+0016. We therefore 
test the sensitivity of the lensing and X-ray analysis to 
the choice of center, finding that the systematic error is 
small. For example, if we adopt the BCG as the cen- 
ter of the X-ray analysis of A1914 then the hydrostatic 
mass estimates change by 1 %, 10%, and 3% at A =500, 
1000, and 2500, respectively. Similarly, if we adopt the 
X-ray center as the center of the lensing analysis of the 
same cluster, then the lensing mass estimates change by 
0.1%, 0.2%, and 0.4% also at A =500, 1000, and 2500, 
respectively. The X-ray analysis of the two clusters was 
therefore revised using the weak-lcnsing determined clus- 
ter center. 

2.4. X-ray morphology 

The X-ray morphology of each cluster was determined 
by Okabe et al. (2010b) by calculating asymmetry (A) 
and fluctuation parameters (F; Conselice 2003) from the 
XMM-Newton X-ray images. In summary, the asymme- 
try parameter is defined BsA = (^2ij\ lij ij ^ i 
the normalized sum of the absolute value of the flux 
residuals. Jjj is the element of the XMM-Newton 
MOS1+MOS2 image in the 0.7-2.0 kcV band, which is 
flat fielded, point source subtracted and re-filled assum- 
ing a Poisson distribution, and binned by 4" x 4". Rij is 
the element of the image derived by rotating the above 
image by 180°. We take into account the position res- 
olution of XMM-Newton by allowing the cluster center 
falling into any neighboring pixel of the r < 4" circle cen- 
tered at the BCG and include this error in quadrature 
in calculating A. Dynamically immature clusters often 
show both an asymmetric X-ray morphology and an off- 
set between weak-lensing and X-ray centers. Therefore A 
is very sensitive to cluster dynamical state. The fluctua- 
tion parameter is given by F = hj - B t3 )/ £V I i3 , 
in which Bij is the element in the smoothed image. This 
parameter describes the degree of deviations from the 
smoothed distribution. We measure the errors of A and 
F assuming a Poisson noise computed within a radius of 
r^QQ , excluding CCD gaps and bad pixels. 

The F versus A plane is divided into four quadrants 
with cuts at A = 1.1 and F = 0.05 in Figure 1 in Ok- 
abe et al. (2010b). Our sample of 12 clusters occupy 
these quadrants as follows: (1) RXCJ2129, A209, A383, 
A1835, and A2390 have both low A and low F; (2) A2261 
and A1914 have high asymmetry parameters; (3) A68, 
RXCJ2337, A267, and Z7160 have high fluctuation pa- 
rameters; and (4) A115 (south) has both high A and 
high F. The clusters falling into quadrant (1) are defined 
as dynamically undisturbed clusters, and the remaining 
clusters as disturbed clusters because high A and/or high 
F indicates that the cluster is still dynamically young. 
It is worth noting that the five undisturbed clusters have 
< 0.06r25oo offset between the X-ray cluster center and 
the weak-lensing cluster center which is the BCG posi- 
tion. 

3. RESULTS 

3.1. X-ray hydrostatic mass versus weak-lensing mass 



The comparison between X-ray and weak-lcnsing 
masses for individual clusters is shown in the left panel 
of Figure [TJ The X-ray to weak-lensing mass ratios vary 
in the range of ~ 0.55 — 1.72. Undisturbed clusters have 
X-ray to weak-lcnsing mass ratios that generally increase 
toward smaller cluster-centric radii (higher over-density). 
In contrast, disturbed clusters show a much greater di- 
versity, with some disturbed clusters having mass ratios 
that increase toward larger cluster-centric radius. The 
distribution of (M x - M WL )/A/ WL at A = 500 (right 
panel of Figure [2]) reveals that the well-known merging 
cluster A1914 is a ~ 5a outlier based on a naive calcula- 
tion of the mean mass ratio for the other 11 clusters. As 
pointed out in Section [2~2l A1914 also shows the largest 
offset (0.44r25oo) between the X-ray and weak-lcnsing de- 
termined cluster centers. 

The average X-ray to weak-lensing mass ratios for 
the full sample and the undisturbed and disturbed sub- 
samples were calculated taking into account the errors in 
the X-ray and weak-lensing mass estimates — see Table[3] 
and the right panel of Figure [TJ M X /M WL is consistent 
with unity across the full radial range probed by the data 
for the full sample of 12 clusters. This also holds for the 
seven disturbed clusters, albeit with uncertainties ~ 2x 
those for the full sample (as seen in simulations; e.g., Na- 
gai et al. 2007). In contrast, undisturbed clusters show a 
gentle decline in M X /M WL to larger cluster-centric radii 
- at A = 500, the five undisturbed clusters show an 
average X-ray hydrostatic mass 9% ± 6% lower than the 
average weak-lensing mass. 

Interestingly, given the apparent trend of M X /M WL 
with radius, Af x /M WL is consistent with unity for the 
full sample at A = 500. As noted above, A1914 — the 
most extreme of the disturbed clusters has a mass ratio of 
- 1.7 at A = 500 (see Tableland Figure [TJ. We there- 
fore investigate the extent to which A1914 might be dom- 
inating the results. We recalculated M X /M WL for the 
full sample and the disturbed clusters excluding A1914 
(right panel of Figure [TJ and Tabic [3]) , finding that the 
average X-ray hydrostatic mass is now lower than the av- 
erage weak-lensing mass just by 6%±5%. Within the un- 
certainties, our result that X-ray hydrostatic and weak- 
lcnsing masses agree for the full sample across the full 
radial range probed, is therefore insensitive to the inclu- 
sion/exclusion of A1914. However, there is evidence for 
shock heating of the ICM in the entropy map of A1914. 
In addition to non-thermal pressure support (likely the 
main reason for X-ray hydrostatic masses being under- 
estimated for undisturbed clusters), there may also be a 
competing effect associated with adiabatic compression 
and/or shock heating of the intracluster gas which leads 
to an overestimate of X-ray hydrostatic masses for some 
disturbed clusters. We therefore argue that robust con- 
straints on the bias of X-ray hydrostatic mass estimates 
for precision cluster cosmology requires both statistically 
large and complete (unbiased) samples in order to sam- 
ple the full range of physical processes at play within the 
underlying cluster population. 

We also calculated the cumulative probability distribu- 
tion function of the X-ray to weak-lensing mass ratio as- 
suming each data point to be a Gaussian distributed vari- 
able and taking into account the errors with 500 Monte 
Carlo simulations. The mean and its standard error are 
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listed in Table HJ Again, a clear trend is found that the 
average hydrostatic to weak-lensing mass ratio declines 
with cluster-centric radius for undisturbed clusters. The 
mass ratios for the full sample and disturbed clusters are 
consistent with unity across the full radial range. 

To further test our results on undisturbed and dis- 
turbed clusters based on our small sample, we applied 
the jackknife method to recalculate the average X-ray 
to weak-lensing mass ratios at A = 500. We randomly 
removed one system from the five undisturbed clusters 
and obtained average X-ray to weak-lensing mass ratios 
in the range of 0.875-0.962. Therefore, our finding that 
the X-ray hydrostatic mass is on average lower than the 
weak-lensing mass for undisturbed clusters is robust. We 
applied the same procedure to the disturbed sub-sample 
and found that the average X-ray to weak-lensing mass 
ratios are in the range of 0.966-1.145. Given the large 
scatter and measurement errors, it is therefore unclear 
whether the average X-ray hydrostatic mass is lower or 
higher than the average weak-lensing mass for disturbed 
systems. 

Despite the modest statistical significance, the anal- 
ysis described above all points toward a trend of im- 
proving agreement between M x and M WL as a func- 
tion of increasing over-density. We therefore proceed 
a simple fit to the data, and obtain the following re- 
lation: M X /M WL = (0.908 ± 0.004) + (0.187 ± 0.010) • 
log 10 (A/500). These results are consistent with those 
of Mahdavi et al. (2008), who found an X-ray to weak- 
lensing mass ratios of 1.06, 0.96, and 0.85 at A = 2500, 
1000, and 500 respectively, with a typical error bar of 
10%. Our trend is slightly shallower than Mahdavi ct 
al.'s result, but is in agreement within the uncertainties. 

Finally, we investigate the issue of scatter. The undis- 
turbed clusters present a factor of ~ 2 less scatter around 
their average X-ray to weak-lensing mass ratio than dis- 
turbed clusters (Table |3j. A key question is whether 
this lower scatter implies lower intrinsic variance within 
the undisturbed sub-sample. We therefore calculated the 
real variance following Appendix A of Sanderson & Pon- 
man (2010) and found that the real variance for disturbed 
clusters is ~ 5 — 10 x larger than that for undisturbed 
clusters at A = 2500, 1000, and 500 (Table 0J). This 
confirms that the smaller scatter measured for undis- 
turbed clusters reflects low intrinsic variance in this clus- 
ter population. This result is in agreement with studies 
based on numerical simulations, and is fully expected 
in both simulations and observations because of the non- 
smoothness of the gas distribution, complex thermal- and 
non-thermal-structure, deviations from spherical symme- 
try, etc., that are typical of dynamically immature clus- 
ters (e.g., Poole et al. 2006; Fabian et al. 2008; Zhang 
et al. 2009). We also stress that the fact that disturbed 
clusters are also typically less spherical than undisturbed 
clusters will also render the deprojection of the lensing 
signal via fitting a three-dimensional NFW model less 
valid in disturbed clusters than in undisturbed clusters. 
A thorough investigation of these issues requires a large 
complete sample of clusters with deep X-ray and lensing 
data. 

3.2. X-ray to weak-lensing mass ratio versus 
morphology indicators 



We now investigate the dependence of the X-ray to 
weak-lensing mass ratio on the morphological parame- 
ters, A (asymmetry) and F (fluctuation), and the con- 
centration of the best-fit NFW halos from Okabe ct al. 
(2010a, 2010b). As the mass estimates for individual 
clusters have large uncertainties, a detailed quantitative 
study of the relationship between X-ray to weak-lensing 
mass ratio and A and F is beyond the scope of that pos- 
sible with the current sample. We therefore restrict our 
attention in this section to general trends that will be 
worth following up with future larger samples. 

First, we plot the X-ray to weak-lensing mass ra- 
tio versus the asymmetry and fluctuation parameters 
in Figure [3J There is a general trend of decreasing 
M X /M WL with increasing asymmetry in the sense that 
the most asymmetric clusters have the largest mass dis- 
crepancies with M x < A/ WL . A1914 is an obvious 
outlier from this apparent anti-correlation. We there- 
fore exclude A1914 and ignore the observational error 
bars in order to fit a "toy-model" to the data, obtain- 
ing M X /M WL = (1.408 ± 0.257) - (0.454 ± 0.263) • A. 
We also note that at fixed asymmetry undisturbed clus- 
ters generally have lower M X /M WL than disturbed clus- 
ters, with undisturbed systems possibly tracing a steeper 
and tighter trend than the latter. However, we caution 
that these results are preliminary — for example, the 
straight-line fit discussed above is dominated by the left- 
most point in the left panel of Figure [31 namely A68. No 
trend is suggested between X-ray to weak-lensing mass 
ratio and fluctuation parameter in the middle panel of 
Figure [3J 

The mass and concentration (c) of an individual NFW 
halo are anti-correlated. For example, the covariance 
for M^qq and c 500 is o 2 M = 1.25, o~Mc = —0.44, and 
a 2 = 1.20 for the sample of 12 clusters. In practical 
terms, this means that clusters with higher wcak-lcnsing- 
based values of C5oo/(c5oo) will on average have higher 
values of M X /M WL , in which (C500) is the average of the 
concentration parameters from Okabe et al. (2010a) for 
the sample of 12 clusters, weighted by the inverse square 
of the errors. Most significantly, the uncertainties on 
c 50o/(c5oo) and M X /M WL are correlated, which may in- 
duce a correlation between the quantities themselves. We 
therefore plot M X /M WL versus C5oo/(c5oo) in the right 
panel of Figure [3J — no obvious relation is seen in this fig- 
ure, suggesting that the effect is small and/or the effect is 
canceled out by other effects. This result does not change 
when we use the concentration parameter c v j r = r v - lr /r s . 
A much larger sample and more detailed analysis are re- 
quired to investigate this issue further. 

Results from numerical simulations indicate that the 
cluster concentration parameter may be related to the 
cluster dynamical state (e.g., Neto et al. 2007). Using 
the Millennium Simulation, Neto et al. pointed out that 
the concentrations of out-of-equilibrium halos tend to be 
lower and have more scatter compared to their equilib- 
rium counterparts. This can also be investigated in the 
observed M X /M WL versus 0500/(0500) plane, however, 
the current sample is too small. 

Finally, we note that the concentration measurements 
used in this article are based on weak-lensing constraints 
(Okabe et al. 2010a). Future papers in this series will use 
more precise concentration parameter estimates based on 



6 



Zhang, Okabe, Finoguenov, et al. 



combined strong- and weak-lensing models of the clus- 
ters. This will allow a more quantitative comparison be- 
tween observations and simulations. 

3.3. Gas mass fraction 

We define the gas mass fraction as Af gas (< 
r^ L )/A/ WL (< r^ L ), and show the gas mass fractions 
for individual clusters at A = 2500, 1000, and 500 in 
the left panel of Figure [4] There is a clear trend of in- 
creasing gas mass fraction toward larger cluster-centric 
radius. A similar trend is also common in both adiabatic 
simulations and in simulations with cooling (e.g., Evrard 
1990, Lewis et al. 2000, Kravtsov et al. 2005). We also 
calculate the average gas mass fractions for the full sam- 
ple and the undisturbed and disturbed sub-samples, at 
A = 2500, 1000, and 500 — right panel of Figure Hand 
Table [5] The average gas mass fraction shows no depen- 
dence on cluster morphology, and approaches the cosmic 
mean baryon fraction f2b/^m = 0.164 ± 0.007 (Komatsu 
et al. 2009) at large cluster-centric radius. At A ~ 500, 
the average gas mass fractions stand at 90% of the cosmic 
mean value, in good agreement with simulations (e.g., 
Nagai et al. 2007). 

Our gas mass fractions are also in good agreement 
with Umetsu et al.'s (2009) joint Subaru weak-lensing 
and AMiBA Sunyaev-Zel'dovich effect (SZE) analysis 
of four clusters. Umetsu et al. obtained (/ gas 2500) = 
0.105±0.015±0.012and (/ ga s,50o) = 0.126±0.019±0.016, 
where the first error is the statistical error, and the 
second one is the standard error from the average due 
to cluster-to-cluster variance. However, the gas mass 
fraction obtained by Mahdavi et al. (2008) using weak- 
lensing and X-ray data is slightly higher than our results: 
(/gas, 2500} = 0.119 ± 0.006. A detailed understanding of 
this ~ 2cr discrepancy awaits analysis of our complete 
volume-limited sample in a future article. 

Several X-ray-only studies have advocated gas mass 
fractions as a potential probe of the dark energy equa- 
tion of state parameter w (Vikhlinin et al. 2006; Allen 
et al. 2008). We therefore investigate how our weak- 
lensing-based gas mass fractions compare with X-ray- 
only gas mass fractions, defining the latter as M gas (< 
r^)/M x (< r^) — i.e., the X-ray-only gas mass frac- 
tions are measured at radii determined from the X-ray 
data derived hydrostatic mass distribution. We cal- 
culated the average X-ray-only gas mass fractions at 
A = 2500, 1000, and 500 for the full sample and the 
undisturbed/disturbed sub-samples (Table[6]). We found 
that the X-ray-only gas mass fractions are statistically 
indistinguishable from weak-lcnsing-bascd gas mass frac- 
tions at all over-densities considered and for all three 
samples, except for disturbed clusters at A = 500, for 
which the X-ray-only value appears to be biased low. 
The most significant result in the context of cluster cos- 
mology is the very close agreement at A = 2500 between 
the respective gas mass fraction measurements for undis- 
turbed and disturbed clusters, regardless of whether the 
cluster mass measurement is based on X-ray or lensing 
data. 

We examine this apparent morphological independence 
of gas mass fraction measurements by looking at the pos- 
sible dependence of gas mass fraction measurements of 
individual clusters on the asymmetry and fluctuation pa- 
rameters. We plot gas mass fraction versus asymmetry 



and fluctuation parameters in Figure [5j The former re- 
sembles a scatter plot as might be expected based on 
the results above. However, undisturbed and disturbed 
clusters both follow the same trend (right panel of Fig- 
ure [5]) of increasing gas mass fraction with fluctuation 
parameter. This may indicate a direct connection be- 
tween the gas mass fraction and substructure fraction, 
i.e., the fraction of cluster mass that resides in mas- 
sive sub-halos (substructures) within the cluster halo. 
Large substructure fractions are generally attributed to 
recent infall of galaxy groups/clusters into more mas- 
sive clusters — see Richard et al. (2010) for more de- 
tails. In a future article, we will investigate the rela- 
tionship between weak-lensing-based substructure frac- 
tions and cluster X-ray morphology, including the fluc- 
tuation parameter. For now, we fit a straight line to all 
12 data points in the right panel of Figure [5l obtaining: 
Ms as /M WL = (0.112 ± 0.011) + (0.502 ± 0.197) • F at 
A = 500. 

Finally, we compare the gas mass fraction with the 
weak-lensing mass. In this current narrow mass range, 
there is no obvious dependence on the weak-lensing mass. 
However, X-ray studies for samples with broad mass 
ranges have shown that there is a strong mass depen- 
dence (e.g., Pratt et al. 2009). This will be discussed 
later in Section |U 

4. DISCUSSION 

In this section, we compare our results with recent nu- 
merical simulations of clusters. There is a broad consen- 
sus among numerical studies that total cluster masses 
derived from X-ray data underestimate the true cluster 
masses in both undisturbed and disturbed clusters for 
variety of reasons, most notably because of the deviation 
from H.E. and extra pressure support (e.g., turbulence) 
beside the thermal gas (e.g., Evrard 1990; Lewis et al. 
2000; Miniati et al. 2001; Miniati 2003; Rasia et al. 2006; 
Nagai et al. 2007; Pfrommer et al. 2007; PV08; Jeltema 
et al. 2008; Lau et al. 2009; Meneghetti et al. 2010). 
An obvious difference between numerical and observa- 
tional studies is that the true cluster mass is known in 
the former, and not in the latter. One solution is to con- 
struct fake weak-lensing data from the simulated cluster 
data (e.g., Meneghetti et al. 2008, 2010). However, large 
samples of such "observed" theoretical clusters are not 
yet available. We therefore rely on the insensitivity of 
the weak-lensing mass measurements to cluster thermo- 
dynamics to assume that our weak-lensing masses are, 
on average, well matched to "true" masses of the nu- 
merical studies. We therefore treat weak-lensing-based 
mass measurements as the "truth" in the comparisons 
described in this section. 

4.1. Summary of recent simulated cluster samples 

Before describing the comparison with simulations in 
the next section, we first summarize the key features of 
the simulated cluster samples against which we compare 
our results. We discuss in turn the samples of Nagai et 
al. (2007), Lau et al. (2009), PV08, and Jeltema et al. 
(2008). 

Nagai et al. (2007) studied 16 simulated clusters with 
T > 2 keV in the mass range of M 500 ~ (0.3 — 9) x 
10 14 M Q /i _1 . They derived hydrostatic masses by ana- 
lyzing mock Chandra data (three projections per cluster) 
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following the method in Vikhlinin et al. (2006), which 
is similar to our reduction method (Zhang et al. 2008). 
We adopt Nagai et al.'s morphological classifications — 
their 48 simulated X-ray cluster observations comprise 
21 undisturbed ("relaxed" in their terminology) clusters 
and 27 disturbed ( "unrelaxed" ) clusters. 

Lau et al. (2009) studied the same simulated clusters 
as Nagai et al., but estimated hydrostatic masses directly 
using the three-dimensional gas profiles. Importantly, 
Lau et al. used mass weighted temperatures (Tmw), in 
contrast to Nagai et al.'s temperatures that were recon- 
structed from mock observations. Lau et al. suggest that 
this difference explains the different average mass biases 
reported in the two papers. The reconstructed temper- 
ature profiles used by Nagai et al. are not spectroscopic 
ones, because they are derived by weighting the contribu- 
tion of temperature components along the line of sight to 
correct for the spectroscopic bias. Nevertheless, they are 
derived from spectroscopic data and sensitive to the local 
multiphase structure of the gas. In contrast to PVge2keV 
(see below) and Nagai et al., Lau et al. define a cluster 
as relaxed/undisturbed only if it appears so in all three 
projections, which gives 6 relaxed/undisturbed clusters 
and 10 unrelaxcd/disturbcd clusters in their sample. The 
different fractions of undisturbed and disturbed clusters 
found by Nagai et al. and Lau et al. (based on the same 
simulated sample) arise from the fact that the perceived 
dynamical state can depend on viewing angle when using 
simple diagnostics. 

PV08 analyzed 100 simulated clusters (three pro- 
jections per cluster) with T > 2 kcV (8.2 x 
IO^MqIi- 1 < M 200 £ 1-2 x W^MqH- 1 ) — this is to 
date the largest, complete volume- limited sample of sim- 
ulated clusters for which the hydrostatic mass biases have 
been investigated in detail. They applied different tech- 
niques to measure the X-ray hydrostatic mass in order to 
disentangle biases of different origin. Here, we compare 
our results with their average biases derived adopting an 
extended /3-model to fit the gas density radial profile to 
estimate the hydrostatic mass, which is similar to the 
procedure in Vikhlinin et al. (2006), Nagai et al. (2007), 
and Zhang et al. (2008), and three-dimensional mass- 
weighted temperature profiles (Tmw) or spectroscopic- 
like temperature (Tsl; Mazzotta et al 2004). 

In order to improve the coverage at high mass end, 
we use an extended version of the PV08 sample, which 
is extracted from a larger cosmological simulation. The 
sample is constructed and analyzed exactly as the PV08 
sample and yields results fully consistent with the lat- 
ter. The new sample (PVgc2kcV sample, hereafter) com- 
prises ~ 120 clusters (3 projections per cluster) with 
T > 2 kcV. Average mass biases derived for the Tmw case 
are fully consistent with those derived using directly the 
three-dimensional gas profiles. The PVge2keV sample 
of 360 simulated X-ray cluster observations were divided 
into 180 undisturbed clusters and 180 disturbed clusters 
using their mock X-ray data; this matches the roughly 
50-50 split between disturbed and undisturbed clusters 
in the observed sample. 

Jeltema et al. (2008) analyzed 61 simulated clusters 
giving a 16% bias of the hydrostatic mass estimate on 
average at A = 500. However, their analysis is restricted 
to just A = 500, and their X-ray hydrostatic masses are 



derived using the true three-dimensional gas profiles only. 
We therefore chose not to include a detailed comparison 
with their results. Nevertheless, notice that their results 
are fully consistent with those presented in this paper. 

In summary, the different hydrostatic mass reconstruc- 
tion methods, and in particular temperature definitions, 
adopted in the simulations discussed above imply that it 
is sensible to compare our observational results with (1) 
PVge2keV (T M w case) and Lau et al. (2009), and (2) 
PVge2keV (T SL case) and Nagai et al. (2007), with an 
emphasis on the latter to match our use of observed spec- 
tral temperatures. We note that the simulated clusters 
extend to lower temperatures (T ;> 2keV) than the ob- 
served sample (T > 5keV). However, the numerical stud- 
ies have shown that the dependence of hydrostatic mass 
bias on cluster temperature (mass) is negligible. The way 
to define the cluster dynamical state in PV08 and Nagai 
et al. (2007) should be more consistent with the one in 
our observations. However, for our comparison there is 
no point to correct the difference because the errors are 
tiny for the samples in simulations compared to the er- 
rors for the observational sample. Finally, we note that 
the cosmological parameters used in the simulations are 
slightly different from those used in this article; the im- 
pact of these small differences on our results/discussion 
is negligible. 

4.2. Comparing observational mass estimates to 
simulations 

We overplot the X-ray hydrostatic to weak-lcnsing 
mass ratios from the simulations in the right panel of 
Figure [TJ We find that there is a general agreement be- 
tween observations and simulations in which the mass 
ratios are consistent with unity at small radii, e.g., at 
A = 2500, and the agreement deteriorates at larger radii. 
Most strikingly, the simulations agree well with our re- 
sult that the mass ratio for undisturbed clusters declines 
gently to larger radii. We estimate that a firm detection 
of the discrepancy between simulations and observations 
for undisturbed clusters, would require a sample of ~ 60 
clusters. 

For disturbed clusters, the simulations — particularly 
those using the spectroscopic-likc temperatures — show 
a significant discrepancy between the X-ray hydrostatic 
and true masses at large radii, while no discrepancy is de- 
tected observationally between the X-ray hydrostatic and 
weak-lcnsing masses. The hydrostatic mass bias seen in 
simulations, when spectroscopic temperatures are used 
in the derivation of hydrostatic masses, might be ex- 
acerbated by over-cooling. Even when small-scale cold 
clumps of gas are masked in the mock analysis, over- 
cooling can lead to underestimated spectral temperatures 
because of the presence of a diffuse cold gas component. 
This may be the cause of the disagreement between the 
observed sample of 12 clusters and the simulated samples 
at large radii. 

4.3. Comparing observational gas mass fractions with 
simulations 

Within a given density contrast, the gas mass fraction 
tends to increase with the increasing cluster mass (e.g., 
Vikhlinin et al. 2006, Gonzalez et al. 2007, Gastaldello 
et al. 2007, Pratt et al. 2009), with the trend being 



8 



Zhang, Okabe, Finoguenov, et al. 



shallower at large radii. This behavior can be explained 
by the scale dependency introduced by cooling (Bryan 
2000), which regulates the amount of cool gas present in 
clusters. Numerical simulations can also reproduce the 
observed behavior of / gas in clusters (Valdarnini 2003, 
Kravtsov et al. 2005, Ettori et al.2006, Kay et al. 2007), 
and strongly support the cooling model. 

Due to the strong mass dependence of gas mass frac- 
tions seen in clusters, a proper comparison between sam- 
ple averages requires that the two samples approximately 
cover the same mass range. This is accomplished by 
extracting a sub-sample of PVgc2kcV clusters that all 
have masses above the lowest mass of observed clus- 
ters: A/500 > 2.6 • 10 14 M Q . This sub-sample contains 
^Vsub = 68 clusters which are subsequently divided into 
34 undisturbed and 34 disturbed clusters based on the 
cluster dynamical state as explained in PV08. The av- 
erage gas mass fractions for all 68 simulated clusters 
and the undisturbed and disturbed simulated clusters are 
compared to the observed gas mass fractions in Figure 01 
To illustrate the strong mass dependence of gas mass 
fractions, we also show in Figure [4] average gas mass frac- 
tions for a more restrictive sub-sample of 32 PVge2keV 
clusters with M 500 > 5 • 1O 14 M . 

The average gas mass fractions of simulated clusters 
agree with those of observed clusters in the sense that 
they do not show any dependence on cluster morphol- 
ogy. The agreement between simulations and observa- 
tions is most striking at A = 2500, where the observed 
and simulated gas mass fractions agree within the uncer- 
tainties when the mass range of the respective samples is 
matched. Overall there is a good agreement within the 
error bars, taking into account that the trend with mass 
is strong and that the mass function of the two samples 
(from simulations and observations) is not exactly the 
same. 

4.4. X-ray and lensing masses in simulated clusters 

As discussed at the beginning of Section 21 large sam- 
ples of mock weak-lcnsing observations of simulated clus- 
ters are not yet available. The largest sample to date 
is that of Mcneghetti et al. (2010), who studied three 
clusters, with three orthogonal projections per cluster, 
to generate a total set of nine mock observations. Their 
X-ray hydrostatic masses are typically biased low by 5%- 
20% due to the lack of H.E. in their simulated clusters, 
which is in agreement with our results. They also found 
the gas mass to be well reconstructed within the region 
where the X-ray surface brightness profile is extracted. 
Their gas mass measurements are independent of the dy- 
namical state of the cluster, with average deviations of 
1% ± 3% at A = 2500 and 7% ± 4% at 200 < A < 500. 
Although Meneghetti et al.'s sample is small, their re- 
sults agree well with our observational results that the 
gas mass to weak-lensing mass ratios are independent 
of the cluster dynamical state, and supports our finding 
that the radius with A = 2500 appears to be the most 
robust radius at which to use cluster gas mass fractions 
for probing cosmological parameters. 

5. SUMMARY AND OUTLOOK 

We have presented a comparison of X-ray hydro- 
static and weak-lensing mass estimates for 12 clusters 
at z ~ 0.2 for which high-quality XMM-Newton and 



Subaru/Suprimc-CAM data are available within the Lo- 
CuSS. Our main results are as follows: 

• For the full sample, we obtain 1 - M X /M WL = 
0.01 ± 0.07 at an over-density of A = 500. We 
also sub-divided the sample into undisturbed and 
disturbed sub-samples based on quantitative X-ray 
morphologies using asymmetry and fluctuation pa- 
rameters. We obtained 1-M X /M WL = 0.09±0.06 
and —0.06 ±0.12 for the undisturbed and disturbed 
clusters, respectively. 

• The scatter around the average X-ray hydrostatic 
to weak-lensing mass ratio for undisturbed clusters 
is half that of the disturbed clusters due to the 
lower intrinsic variance among the population of 
undisturbed clusters. 

• For disturbed clusters, it is unclear whether the X- 
ray hydrostatic mass is consistent with the weak- 
lensing mass, due to large scatter. In addition to 
non-thermal pressure support, there may be a com- 
peting effect associated with adiabatic compression 
and/or shock heating of the intracluster gas which 
leads to overestimate of X-ray H.E. masses for dis- 
turbed clusters. The most prominent example of 
this in our sample is the famous merging cluster 
A1914. 

• Despite the modest statistical significance of the 
mass discrepancy in the undisturbed clusters, we 
detect a clear trend of improving agreement be- 
tween M x and M WL as a function of increasing 
over-density, M X /M WL = (0.908±0.004) + (0.187± 
0.010) -log 10 (A/500). 

• There is a general agreement between the X-ray 
hydrostatic to weak-lensing mass ratio in observed 
and simulated clusters in which the mass ratios are 
both consistent with unity at small radii (i.e., at 
A = 2500), with the agreement deteriorating at 
larger radii, i.e., out to A = 500. This deterioration 
is dominated by disturbed clusters, with the undis- 
turbed simulated clusters reproducing well the ob- 
served gentle decline in M X /M WL as a function of 
increasing cluster-centric radius. 

• The weak-lensing mass-based cumulative gas mass 
fraction increases with radius, but still lies below 
the cosmic baryon fraction at the largest cluster- 
centric radii probed (i.e., A = 500). An important 
finding is the absence of dependence of the gas mass 
fraction on cluster dynamical state (i.e., X-ray mor- 
phology). The X-ray-only gas mass fractions are 
also consistent with the weak-lensing mass-based 
gas mass fractions at A = 2500, supporting the 
proposal to use this measurement as a probe of the 
dark energy equation of state parameter w. 

In summary, our results demonstrate that XMM- 
Newton and Subaru are a powerful combination for cal- 
ibrating systematic uncertainties in cluster mass mea- 
surements and suggest an encouraging convergence be- 
tween X-ray, lensing, and numerical studies of cluster 
mass. Nevertheless, our observational sample remains 
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very small at just 12 clusters. Our detailed results are 
therefore vulnerable to inclusion/exclusion of extreme 
clusters, as highlighted by our discussion of A1914. Ro- 
bust constraints on systematic uncertainties in cluster 
mass measurement therefore await a thorough investiga- 
tion of a large and complete volume-limited sample of 
clusters, such as that planned within LoCuSS. On the 
theoretical side, it is important to move as rapidly as 
possible toward large samples of mock weak-lensing and 
X-ray observations of simulated clusters, such as those re- 
cently pioneered by Meneghetti et al. (2008, 2010). With 
both of these observational and theoretical data sets in 
hand, a detailed and statistically robust comparison will 
be possible, helping to calibrate future cluster cosmology 
experiments. 
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A A 

Fig. 1. — Left panel : X-ray hydrostatic to weak-lonsing mass ratios as a function of density contrast for individual clusters, in which red 
dashed and blue solid lines denote undisturbed and disturbed clusters, respectively. Right panel: X-ray to weak-lensing mass ratios for all 
clusters (bottom), undisturbed clusters (middle), and disturbed clusters (top) from our sample including (filled black boxes) and excluding 
(open black boxes) A1914, together with the X-ray hydrostatic mass to true mass ratios from numerical simulations (PVge2keV, red circles 
using TmWi blue diamonds using Tsl; Lau et al. 2009, green triangles using quasi- Tmw; Nagai et al. 2007, light blue stars using Tsl). 
The data points at each density contrast are off by 0.022 dex for clarity. We also show the best fit of the X-ray hydrostatic to weak-lensing 
mass ratio as a function of density contract for the five undisturbed clusters, M X /M WL = (0.908 ± 0.004) + (0.187 ± 0.010) ■ log 10 ( A/500). 
(A color version of this figure is available in the online journal.) 
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Fig. 2.— Normalized histograms of (M x - M WL )/M WL for all 12 clusters at A = 2500 (left panel), A = 1000 (middle panel), and 
A = 500 (right panel). 
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Fig. 3. — X-ray hydrostatic to weak-lcnsing mass ratio at A = 500 vs. asymmetry parameter (left panel), fluctuation parameter (middle 
panel) , and C5oo/{c5oo) (right panel) , respectively. The colors and symbols are the same as shown in the left panel of Figure [T] The best 
fit with equal weighting of all data points except for A1914 (the top rightmost point) is shown in the left panel. (A color version of this 
figure is available in the online journal.) 




o 2000 1000 500 2000 1000 500 

A A 

Fig. 4. — Left panel: gas mass fractions using weak-lensing masses for individual clusters. The colors, lines, and symbols are the same 
as shown in the left panel of Figure[T] Flight panel: average gas mass fractions using weak-lensing masses (filled symbols). Boxes, circles, 
and triangles denote all, undisturbed, and disturbed clusters. The data sets are off by 0.048 dex at each density contrast for clarity. The 
gray horizontal band shows the lcr range of the cosmic mean baryon fraction from WMAP five-year data, Q^/Qm = 0.164 ± 0.007, in 
Komatsu et al. (2009), and the black horizontal line corresponds to 0.9 X the mean cosmic baryon fraction. Also shown are the gas mass 
fractions for the sub-samples drawn from the PVge2keV sample with mass cuts of M 500 > 5 x 1O 14 M (small green open symbols) and 
A^500 > 2.6 X 10 14 Mq (big magenta open symbols). (A color version of this figure is available in the online journal.) 
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Fig. 5. — Gas mass to wcak-lcnsing mass ratio at the radius with A = 500 vs. asymmetry parameter (left panel) and fluctuation 
parameter (right panel). The best fit with equal weighting of each data point is shown in the right panel. The colors and symbols are the 
same as shown in the left panel of Figure [T] (A color version of this figure is available in the online journal.) 



TABLE 1 

Cluster centers from weak-lensing and X-ray analysis 



Cluster Name 


X-ray Center [J2000] 


Weak-lensing Center 


[J2000] 


Offset 




X-ray Morphology 




R.A. 


Dccl. 


R.A. 




Decl. 


(arcmin) 


^2500 




A68 


00 37 06.2 


09 09 28.7 


00 37 06.9 


09 


09 24.5 


0.18 


0.10 


Disturbed 


A115 (south) 


00 56 00.3 


26 20 32.5 


00 56 00.3 


26 


20 32.5 


0.00 


0.00 


Disturbed 


A209 


01 31 52.6 


-13 36 35.5 


01 31 52.5 


-13 


36 40.5 


0.08 


0.04 


Undisturbed 


A267 


01 52 42.0 


01 00 41.2 


01 52 41.9 


01 


00 25.7 


0.26 


0.14 


Disturbed 


A383 


02 48 03.3 


-03 31 43.6 


02 48 03.4 


-03 


31 44.7 


0.03 


0.01 


Undisturbed 


A1835 


14 01 01.9 


02 52 35.5 


14 01 02.1 


02 


52 42.8 


0.13 


0.06 


Undisturbed 


A1914 


14 26 00.9 


37 49 38.8 


14 25 56.7 


37 


48 59.2 


1.22 


0.44 


Disturbed 


Z7160 


14 57 15.2 


22 20 31.2 


14 57 15.1 


22 


20 35.3 


0.07 


0.05 


Disturbed 


A2261 


17 22 26.9 


32 07 47.4 


17 22 27.2 


32 


07 57.1 


0.34 


0.13 


Disturbed 


A2390 


21 53 37.1 


17 41 46.4 


21 53 36.8 


17 


41 43.3 


0.08 


0.03 


Undisturbed 


RXCJ2129.6+0005 


21 29 39.8 


00 05 18.5 


21 29 40.0 


00 


05 21.8 


0.07 


0.04 


Undisturbed 


RXCJ2337.6+0016 


23 37 37.8 


00 16 15.5 


23 37 39.7 


00 


16 17.0 


0.48 


0.25 


Disturbed 



Note. Within the sample. A1914 and RXCJ2337. 6+0016 show extreme offsets between the X-ray and weak-lensing determined cluster 
central positions. Therefore, the X-ray analysis of the two clusters was revised using the weak-lensing determined cluster centers. 
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TABLE 2 

Cluster mass measurements' 1 



Cluster Name 




M WL 




-WL 

T 2500 
M X 




Ms as 






m wl 




T .WL 
'1000 

M x 




MS 


as 






A/ WL 




' 500 

M x 








A68 


1 


42±[j 


59 
02 


1 


91 


± 





57 





158 


± 





007 


2 


78l° 


vy 

78 




56 


± 


1 


06 





360 


± 


o 


026 


1 


15 


+ 1 
- 1 


■21 
00 


5 


50 


± 


1 


65 





579 


± 





057 


A115 (south) 


1 


22±° 


or, 
07 





68 


± 





12 





107 


± 





005 


2 


5ii° 


!)3 
91 


1 


98 


± 





50 





329 


± 





017 


3 


85 


f l 
-i 


01 
32 


3 


84 


± 





99 





679 


± 





037 


A209 


2 


18±8 


40 

40 


1 


95 


± 





55 





208 


± 





014 


5 


241° 


74 

72 


1 


41 


± 


1 


40 





557 


± 





050 


8 


81 


+- 1 
-i 


31 
20 


6 


57 


± 


1 


95 





963 


± 





101 


A267 


1 


43±g 


25 
25 


1 


66 


± 





15 





139 


± 





007 


2 


4lt°o 


42 

40 


2 


58 


± 





86 





281 


± 





021 


3 


29 


fO 
-0 


08 
01 


3 


75 


± 


1 


16 





426 


± 





039 


A383 


1 


76±C 


21 
2(1 


1 


61 


± 





48 





127 


± 





008 


2 


641° 


43 
40 


2 


42 


± 





72 





227 


± 





022 


3 


38 


+ 
-0 


71 
02 


3 


23 


± 





95 





340 


± 





041 


A1835 


2 




57 
58 


2 


98 


± 





89 





423 


± 





023 


6 




1)5 
90 


5 


66 


± 


1 


68 





844 


± 





080 


9 


65 


H-l 
-1 


70 
51 


8 


59 


± 


2 


50 


1 


261 


± 





154 


A1914 


2 


09±C 


31 
31 


2 


78 


± 





76 





268 


± 





015 


3 


351°, 


511 
47 


1 


36 


± 


1 


22 





497 


± 





038 


1 


46 


+ 
-0 


75 

09 


7 


69 


± 


2 


24 





754 


± 





066 


Z7160 







38 
42 


1 


17 


± 





35 





129 


± 





004 


1 


75±° 


58 
55 


1 


80 


± 





52 





251 


± 





014 


2 


61 


fO 
-0 


97 
82 


2 


16 


± 





72 





401 


± 





033 


A2261 


3 


55±C 


44 

44 


2 


77 


± 





75 





266 


± 





026 


5 


951°, 


74 
71 


1 


54 


± 


1 


20 





516 


± 





065 


8 


12 


+ 1 
-1 


23 
12 


6 


24 


± 


1 


65 





786 


± 





114 


A2390 


3 


14±B 


44 

43 


3 


92 


± 


1 


15 





395 


± 





031 


5 


221° 


82 
77 


6 


02 


± 


1 


91 





785 


± 





083 


7 


09 


H-l 
-1 


2!) 
17 


7 


48 


± 


2 


22 


1 


180 


± 





150 


RXCJ2129.6 


1 


38±C 


53 
54 


1 


75 


± 





52 





166 


± 





009 


2 


971° 


09 
71 


3 


07 


± 





93 





364 


± 





031 


1 


68 


+ 1 

-0 


(19 
97 


4 


46 


± 


1 


29 





581 


± 





065 


RXCJ2337.6 


2 


42±g 


30 
37 


1 


97 


± 





44 





202 


± 





015 


3 


721° 


5(1 
47 


3 


16 


± 





68 





418 


± 





045 


1 


85 


fO 

-11 


75 
70 


4 


15 


± 





97 





652 


± 





084 



Note. a The values are in units of 10 14 solar mass. 



TABLE 3 

Comparison of X-ray and weak-lensing mass estimates 



Sample 


-WL 
r 2500 


M X /M WL 

.WL 
7 1000 


T 500 


All 


1.01 ±0.07 


0.97 ±0.05 


0.99 ±0.07 


Undisturbed 


1.04 ±0.08 


0.96 ±0.05 


0.91 ±0.06 


Disturbed 


0.98 ±0.12 


0.97 ±0.09 


1.06 ±0.12 


A11-A1914 


0.97 ±0.07 


0.94 ±0.05 


0.94 ±0.05 


Disturbcd-A1914 


0.92 ±0.11 


0.91 ±0.07 


0.97 ±0.08 



TABLE 4 

Results of Monte Carlo simulations' 1 



Sample 




M X /M WL 






-Wl 


-Wl 


-WL 




r 25D0 


r 1000 


r 500 


All 


1.03 ±0.09 ±0.13 


0.94 ±0.06 ±0.05 


0.95 ±0.07 ±0.08 


Undisturbed 


1.02 ±0.05 ±0.02 


0.94 ±0.04 ±0.01 


0.90 ±0.04 ±0.02 


Disturbed 


1.05 ±0.11 ±0.22 


0.94 ± 0.07 ± 0.08 


1.02 ±0.08 ±0.11 



Note. a The valu es qu oted for each comparison are the mean, standard error, and real variance based on 500 Monte Carlo simulations 
described in Section |3. II 



TABLE 5 

Weak-lensing-based gas mass fractions 



Sample 




M gzs/ M WL 






,W1 


-WL 






'2500 


Mono 




All 


0.099 ± 0.008 


0.118 ±0.007 


0.130 ±0.008 


Undisturbed 


0.101 ±0.015 


0.118 ±0.012 


0.124 ±0.011 


Disturbed 


0.096 ± 0.009 


0.119 ±0.009 


0.136 ± 0.011 



TABLE 6 
X-ray-only gas mass fractions 



Sample 


r 25nn 


'"innn 


r 5nn 


All 


0.101 ±0.010 


0.118 ±0.015 


0.131 ± 0.015 


Undisturbed 


0.100 ±0.017 


0.122 ±0.028 


0.124 ±0.023 


Disturbed 


0.101 ±0.018 


0.115 ±0.024 


0.110 ±0.026 



